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Within the framework of the QCD sum rules, we consider the three-point correlation function, 
work at the hmit — s- and 0, and pick out the singular term ~ ^ to extract the pionic 

coupling constants of the 1""*" hybrid meson. Then we calculate the decay widths of different modes. 
The decay width of the 5*- wave modes fcivr, /ivr increases quickly as the hybrid meson mass and decay 
momentum increase. But for the low mass hybrid meson around 1.6 GeV, the P-wave decay mode 
p-K is very important and its width is around 180 MeV, while the widths of rjn and tj'tt are strongly 
suppressed. We suggest the experimental search of 7ri(1600) through the decay chains at BESIII: 
e'^e~ J/ijj{tp') ^ TTi -I- 7 or e^e~ J/'ip{tp') —^ttj+p where the tti state can be reconstructed 
^ ■ through the decay modes tti ^ p7r — > ■k'^tt'tv'' or tti — 5- /i (1285)71". It is also interesting to look for 

TTi using the available BELLE/BAB AR data through the process e^e^ — > 7* pyri, fciTri, 77ri etc. 
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I ■ I. INTRODUCTION 

«D : . . nn . 

, A hybrid meson state contains one valence quark, one valence anti-quark and one valence gluon [IHJ], which 

is in contrast with the qq meson in the conventional quark model. Especially the hybrid states with J^'^ = 
0^^, 0^ , 1 ^, 2^ , • • • are of particular interests since these exotic quantum numbers arise from the manifest gluon 
I degree of freedom and can not be accessed by the qq meson. Experimental confirmation of the exotic hybrid mesons 
' and glueballs will be a direct test of quantum chromodynamics in the low energy sector. In the past three decades 
. there have been a lot of experimental and theoretical investigations of the hybrid states. However, their nature remains 
' elusive. 

^ ! Up to now there are three candidates with J^<^ = 1"+: 7ri(1400), 7ri(1600) and 7ri(2000). According to PDG, 
d ■ their masses and widths are (1376 ± 17, 300 ± 40) MeV, {1653tll, '^'^^ttl) MeV and (2014 ± 20 ± 16, 230 ± 21 ± 73) 
' MeV Q. Many collaborations reported evidence of the 1 states such as the CLAS Collaboration, the Crystal 
Barrel Collaboration, the COMPASS Collaboration, the E862 Collaboration etc.. For example, the 7ri(1400) state 
was observed in the reactions n~p — >■ rjTr^n — > n'^TT^r], pn — tt'tt^t] Q, it~p — >■ r]TT~p Q. The 7ri(1600) state 
was observed in the reactions n^p — >■ r]'TT~p tt^jj — ^ ujtt~'k^p fioj and n~p — >■ ?77r+7r~7r~p [ll|. The 7ri(2015) 
meson was observed in the channels n^p — > tJ7r~7r'^p fiol and n^p — > rjiT^TT^TT^p [ll|. Some other recent experimental 
papers include Refs. [l^ - [T5j . We note that the Ref. |15| studied the 7r~7r^7r+ final state whose data show a significant 
' production of 7ri(1600) decaying to pw. We also note that the existence of 7ri(2015) is still questionable, and it is not 
included in the most recent PDG 2010 [l6| . 

There have been many theoretical calculations of the hybrid meson mass in literature p7l - l23l. The 1 ^ mass 
extracted from the quenched lattice QCD simulation ranges from 1.74 GeV [11] and 1.8 GeV (29l| to 2 GeV [27|. 
Among the various phenomenological models, the flux tube model is a popular one. Within this framework the 
glue degree of freedom is modeled by a semi-classical color flux tube. The hybrid meson with the quantum numbers 
JP^ - 1- + was found to be around 1.9 GeV j25l.[26j. The mass from constituent gluon model is around 1.93 GeV [3Q| . 
Many authors studied the mass of the hybrid meson using the QCD sum rule formalism In our previous 

papers, we used the tetraquark currents (qqqq) to study the 1 ^ mesons using the method of QCD sum rule since 
both the tetraquark and hybrid interpolating currents may couple to the same states. The extracted mass is around 
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1.6 GeV and 2 GeV, quite close to the 7ri(1600) and 7ri(2000) ^,\^. 

To learn more about the hybrid states, it is equally important to study their decay properties besides their mass 
spectrum. In the flux tube model, a hybrid meson decays as the flux tube breaks. The "^Pq pair creation" mechanism 
is introduced to create a quark-antiquark pair with S = 1, L = 1 and total angular momentum J = 0. With this 
assumption, the lowest-lying hybrid state with J^*-^ — 1 cannot decay into two ground states, such as tttt, tti] and 
TTp, etc in the original flux tube model [4l|. Later some authors modified this phenomenological model further with the 
introduction of a new decay vertex in order to study the hybrid meson decay process l43j . With this modification, 
the contribution of the P-wave decay mode pn is not negligible. There are also other theoretical approaches on the 
decay, photo- and electro-production of the 1 hybrid meson [35l l43 - l52| . 

In literature, the three-point correlation function was invoked to discuss the decay width of the modes pir, rjn 
etc m, HBj, where the sum rules were derived at the symmetric point = {p' + qY — — —Q^ > 0. Later, the light 
cone QCD sum rules was employed to calculate the decays of the hybrid mesons by one of the present authors [46^, '4^. 

In order to perform a systematical study of the 1 hybrid meson we use a different formalism within the framework 
of QCD sum rule in this paper, which was first developed in Rcf. [sl]. After calculating the three-point correlation 
functions using the method of operator product expansion, we ignore the small pion mass term ra^ in the denominator 
and pick out the divergent term ^ in the limit — > where q is the pion momentum. After comparing the 

singular term of the three-point correlation function both at the phenomenological and quark-giuon level, we make 

Borel transformation of the variables p^,p'^ = {p + qY, and extract the coupling constants (such as gp-n)- Then we 
calculate the decay widths of both the isovector hybrid states tti of I'~^J^'~"' ~ 1^1 ^ and the isoscalar ones cri of 
jGjPC _ ^j^g non-vanishing decay modes include tti — >■ pn, rjTr, ri'ir, bin, fiir, and cti — >■ 7777', citt, /177, etc. 

Our paper is separated into several sections according to the different decay modes. In Sec. [Ill we study the decay 
mode TTi — pn. In Sec. lIIIi we study the decay modes tti — > rjir and t/'tt. In Se. lIVl we study the decay mode tti — > biir 
using the derivative current. In Sec|Vl we study the decay mode tti — )• fin. In Sec lVIi we extend the same formalism 
to study the isoscalar hybrid state with I^J^"^ — 0+1 '". Sec. IVIII is the summary. In Appendix [K[ we study the 
decay mode ni — >■ bin, using the tensor current for bi. 



II. THE DECAY MODE tti ^ pn 
A. Three-Point Correlation Function 

For the past decades QCD sum rule has proven to be a very powerful and successful non-perturbative method [56l . [57l |. 
In the usual QCD sum rule, we consider the two-point correlation function in order to extract the hadron masses. In 
this paper we consider the three-point correlation functions in order to study the decay properties of hadrons: 



The hybrid state tti has several decay modes, such as pn, fi{1285)n, &i(1235)7r etc.. So we need calculate several 
three-point correlation functions. The procedures are more or less the same. In this section we study the decay mode 
TTi pn. First we show the three currents used in our calculation. The hybrid current with J^^ — 1 ^ is: 



^TTl ah ^ /^n h /•o\ 

"97 —9s<^f,i,<l , (^) 



where the summation is taken over repeated indices (/u, v, ■ ■ ■ for Lorentz indices, and a,b, ■ ■ ■ for color indices). The 
hybrid state tti is an isovector state. We study the neutral one, whose quark configurations are {uu — dd)g: 

1 _ A" 

and it couples to 7ri(1600) through 

(0|r;^|7ri(p,A)) =^/2^,r<e^. (4) 
For the p meson, we use the vector current 
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and it couples to the vector meson p through 

(Ob-^|p(p,A)) =mp/pe^. (6) 

For the pseudoscalar meson tt, we use 

jf =dl5U, if = (U75U - ^75^) , il =u-fbd, (7) 

and it couples to vr through 

(0bI|7r(p))=/; = ^. (8) 

Having defined all these interpolating currents and their couplings to the physical hadrons, we can write down the 
three-point correlation function at the phenomenological side: 



gp7Tef,'iy'ai3q p^{gp.p.' 2~)'^9<^iy' j" 

mz, mi 



2 ■ n^2 (mj^ — — ie){m? — p'^ — ie){m'^ — 1^ — ie) 



where the momenta of tti, p and tt are p^, p'^ and q^, respectively. This is for the decay of ttJ — >■ p^7r+, and we can 
obtain the same result for Trf — )• p+Tr^. The coupling constant gp,^ is defined as 

C - gp.e^uap^Td'^^^d^P'" + ■■■ ■ (10) 

To simplify our calculation, we work at the massless pion pole. Such a formalism was first applied to calculate the 
pion nucleon coupling constants very successfully decades ago [ssl - lssl \5T\ . Later the same formalism was employed 
to discuss various pionic couplings [58]. We first take the small tt mass in the denominator to be zero. Then we work 

at the limit 0. T^^^^^'' is divergent at q^ = up to the order of (g^)^^. Therefore, we will also choose such 

divergent terms at the QCD side in the next subsection. 



B. Operator Product Expansion 

In the previous subsection we have obtained the expression of the three-point correlation function at the phe- 
nomenological side. In this subsection we will calculate this at the QCD side using the method of operator product 
expansion (OPE). To do this, first we insert the three currents into the three-point correlation function. After doing 
the Wick contractions we obtain: 

{OlTji;- ix)r\y)^lmO) (11) 



2V2 



l=Tr(^5r(x)7,,'^5^''(-2/)75^^f (2/ - ^h.)Kb9sG;^'{0) . 



2V2' 

Inside this expression, Su and Sd are the quark propagators for up and down quark, respectively. In the presence of 
quark and gluon condensates their expressions are: 

iS'^\x) ^ {0\T[q-{x)q\0)]\Q) (12) 
iS'''' , i X^i, „„ 1 , (5''^_ , S'^^x^ 



g^GZ^-ia^'^x + xaf^n - —{qq) + ——{g^qaGq) , 



27r2a;4 327r2 2"^' ^"'x^' ' 12'^^' 192 
where the masses of up and down quarks have been neglected already. 

The OPE calculation is largely simplified when we work at the pion pole and only choose the terms divergent at 
the q'^ ^ hmit: 

f{gsqaGq),3 , 1. {qq){g^,GG) , 1 , l.X 

iP^P^D- q2 [ V 18^2 V P'^^^ ^ ' 

We note here that the tri-gluon condensate (g^fGGG) as well as the correction vanishes in this case. 



4 



C. Numerical Analysis 

In our numerical analysis, we use the following values for various condensates and mg at 1 GeV and as at 1.7 GeV 

{m) = -(0.240 GeV)^ 
(ss) = -(0.8 ± 0.1) X (0.240 GeV)^ , 
{glGG) = (0.48 ± 0.14) GeV'' , 
{glGGG) = 0.045 GeV^ , 

{gsqcjGq) = -Ml X {qq) , (14) 
Afo^ = (0.8 ± 0.2) GeV^ , 
ms{l GeV) = 125 ± 20 MeV , 
as(l-7GeV) = 0.328 ± 0.03 ± 0.025 . 

There is a minus sign in the definition of the mixed condensate (ggqaGq), which is different from that used in some 
other QCD sum rule studies. This difference just comes from the definition of coupling constant gs ^59j . We use the 
following values for the initial and final states [U [E^, EHl : 

m^, = 1.6GeV , = 0.026GeV , (15) 

™. = i40Mev,/; = ^ = Mz^l^M£X)!, 

U 131MeV 
nip = 770MeV , fp = 220MeV . 

We can compare the three-point correlation function at the phenomenological side Eq. ([9]) and at the QCD side 
Eq. (I13|) in order to calculate the coupling constant gp^^. Due to p = p' -\- q, we have several different choices for the 
QCD sum rules analysis. Here we use p'^ and q^, p^, p'"^ and q^ instead of pp, p ■ p' , p ■ q and p' • q. 

In the QCD sum rule, we use the Borel transformation to suppress the high order contributions. We use T to 
denote the Borel Mass Mb- Here we have two types of Borel transformation: 

1. Assume p^ = p'^ and perform the Borel transformation once, B{p^ = p'^ — > T^); 

2. Perform the Borel transformation twice, B{p^ — ^ T^, p'^ — ^ T2), and then assume Ti — T2 — T. 




I . . . . ^0 

0.5 1 1.5 2 2.5 3 
Borel Mass- [GeV^] 

FIG. 1: The coupling constant Qp^ as a function of Mg, calculated using Eq. (|16p . 



For the decay mode pn, we employ the first type of Borel transformation. We compare Eqs. ^ and (jl3p . assume 
p2 _ p/2 perform the Borel transformation once, B{p^ = p'^ T^) to arrive at 

V^U.mljpmpfi^ ( ,„^2 „„2/j,2\ _ 2{gqaGq) {qq){glGG) 1 

ml-ml^ r e . J - . (lb) 

Using this equation, the coupling constant can be calculated as a function of Af^, as shown in Fig.[T] The result 
is around 4 GeV^^. The formula of the decay width reads 

r(7r?^p+^-+p-7r+) = 2x||l|g;|3, (17) 
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where (f^r is the momentum of the final state tt: 



[(7774 -(mp + m^)2)(m2^-(mp-m^)2)]i/2 
The decay width of tti — >■ pir is around 180 MeV. 

III. THE DECAY MODES vn rjTv,r]'-K 

A. Three-Point Correlation Function 

In this section we study the decay modes tti rjir and tti — > rj'ir. The interpolating current for the 77 meson is 

j'' = cos6'pj''«-sin0pj''°, (19) 
where the mixing angle is 9p = —17° [69j . and it couples to r/ through 

(0|/'h(p))=A,. (20) 
The S'f/(3) octet and singlet currents in Eq. (IT51) are 

j''^ = -i=(u75U + d75rf-2s75s), (21) 

j'"' = ^(u75U + ^75*^ + S75S) . 
v3 

The interpolating current for the rj' meson is 

f = sin Opf'' + cos Opf , (22) 

and it couples to rj' through 

(0|/h'(p))=A„. (23) 
The three-point correlation function for the decay mode tti rji: is: 

T;"^™Hp,p',g) - y'dWye*^''^e*««(0|Tj''(a:)j-(y)77t(0)|0) (24) 

I3^.9a + g^^PaAffM- J (^2^ _ p2)(^2 _ p,2)(^2 _ ^2) ' 

= -(o^n +0^11' ^ A PmP • g „B P,^P ■ P' ^ V^/^i ^^^l^\f'. 

(9^.1^ + y,,./'^ ^2^ ff^- ^2^ ^ (^2^ _ p2)(^2 _ p/2)(^2 _ ^2 ) ' 

where the coupling constants g^^ and g^^ are defined in the following Lagrangian 

^ = W^7rT^iidt^'^)v + ig^Tr'^iidf^v)'^ ■ (25) 

In order to calculate the three-point correlation function tJJ^^^^"^ at the QCD side, once again we work at the pion 
pole and only choose the terms divergent at the — > limit at the quark and gluon level. When — > 0, we obtain 



q' 



24V27T^ Jo 2(ri+T2) 
{glfGGG) f°° 2^ ,2 / 2ti , -t2 



■■T;!<^^^\p,p',q) ^ X{e) X (,^,(p2_y2)M^ re-^^W^--^dndr2 (26) 



{gJGGGl / r^p^+r2p'^ ( 2Ti 



96V27r2 7o V^^(ri-|-r2)2 ' ^^(ri+r2)2 

(5?/GGG) 



1/2 I2\ ('''l + '''2 )'''2 ,1(2 a\ ''"2 A J J 



48%/27r2 y ' 
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where X(9) — cos9p ^ y | sin^p. In the calculations, we have used the following equation [67 



(a; - y)2'2/2m^2« A'+"'^+'^-illmln\ J„ (riT2 + r2r3 + T3ri)'+™+"-2 
We find that only when i ^ j, the term 



(27) 



Iq iriT2 + T2T3 + T^Tiy 

is divergent at (7^ — > (it is up to the order of (q^)-'^^^^)- 



B. Numerical Analysis 

To perform the numerical analysis, we use the following values (65j 

m„ = 0.547GeV , A„ = 0.23GeV2 . 



(28) 



0.958GeV,A,,/ = 0.33GeV^ 



There are two independent Lorentz structures, and p'^. We will use both of them to perform the QCD sum rule 
analysis. 



1. Lorentz structure 

First we choose the Lorentz structure q^. If we perform the Borel transformation once, the condensate {g'^GG) 
vanishes, and only the condensate {g^fGGG) survives. Consequently, the coupling constants is tiny. 
If we perform the Borel transformation twice, the condensate {g^GG) does not vanish, and we obtain 



(29) 



X{9){glGG) ( _d d T2 \ 



24V27r2 V ^dn dT2'2{Tl+T2Y 

X{e){glfGGG) ( 2ti ,d d Arl + t^)t2 



96V27r2 Mn+ra) 



-dri dT2' (ri + T2) 



ri=l/Tf,r2 = l/T= 



Using Eq. ([29)) we can calculate [g^^ — g^^). The coupling constant {g^,^ — g^'j,) can be obtained similarly. They 

1 2 




1 1.5 2 2.5 
Borel Mass^ [GeV^] 



0^ 
0.5 



1 1.5 2 2.5 
Borel Mass- [GeV^] 



1.5 
1 

0.5 




FIG. 2; The coupling constant gr^-w and g^/^ as functions of Mg. The left figure is for {g^^ — (?^^), calculated using Eq. ([29] 



and the right figure is for (g^/^ — g^/.,^). The Lorentz structure here is g^. 



are functions of Mb, as shown in Fig. [5J The L.H.S. is for {g^^ — g^^^), and the result is around -0.7; the R.H.S is for 
((7^,^ — g^iyr), and the result is around -0.4. Using the formula of the decay width 



F« ^ 7/^0) 



247r ml 



(30) 



the decay width of tti r/7r is around 0.9 MeV, and the one of tti — >■ t^'tt is about 0.1 MeV. 



7 



2. Lorentz structure 

Next we move on to the Lorentz structure p'^. After performing the Borel transformation twice, we obtain 



X{e){glfGGG) / -r2 ,d d . t| 

96y27r2 V (ti + Ta)^ ^ Oti dT2 ' (n + T2, 



ri=l/T2,r2 = l/T| 



However, the condensate (g^GG) vanishes here, and the obtained resuh is very smaU (around -0.05 for (g^^ — g^^)). 

IV. THE DECAY MODE TTi ^ 6i(1235)7r 
A. Three-Point Correlation Function 

In this section we study the decay mode tti — > foi(1235)7r. For the fei(1235) meson, we use the following current 

jl' = {dDl-f^u) = --i=(ui?^75M - d^la^d) = (ul?^75d) , (32) 
and it couples to 6i(1235) through [E^l 

(0|j^M&i(1235)(p,A)) = /,,e^. (33) 
Then we can write down the three-point correlation function at the phenomenological side: 

r^'r*™Hp,p',g) - J d^xd^ye^P'^e^'^y{0\Tjt^{x)r\y)vlmO) (34) 

mi^ m^^ [mi^ - p^)[mf^ - p'^}{mi - q^) 

X {9H^9^i'>y' + 91 {{p ■ p')g^i'u' ~ p'f,'Pu') + 92 {{v ■ q)9t,'y' - (l^i'Py') 

+9f{{p' ■ (i)9t^'y' -p'^'qw) + gfiq'^gf.'u' - v?!.')) • 

This is for the decay mode ttj* b^n'^. The same result holds for tt" b^ir^ . The coupling constants g^_^^ and gf 
(1=1,2,3,4) are defined through the following effective Lagrangians: 

£s = 9lA X hp ■ ^9"^^ , (35) 

>Cd = +.9?(5aTfl^ - a^Tfl„) X d'^bl ■ TT (36) 
+92{daT^ip - dfiTTia) X b1 ■ (9"7f 

+.9?7ff X (9a6i/3 - dpbia) ■ d^n 

+gfiTTif3 X 6f • a„a"7r - 7?!^ X ■ d'^d^n) , 

Once again, the OPE calculation is largely simplified when we work at the pion pole and only choose those divergent 
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terms at — ^ 0. When 0, we obtain 

2 



.2 

2 



, 6riT2 + 4r| , -4t2 , , -2r| 



, / 2 a\ ^2 , I 2 /2^ (~2ri + 2rir2 - 2t|)tiT2 
.11(2 /2N -2Tf \ {glfGGG) Pf,p^ 

+PuP^ [P ~P ) 7 ^ 73 '^ndT2 ^ „ ^ . 

We note that we do not inchide the higher order as corrections such as as{gsGG), etc. 



B. Numerical Analysis 



For numerical analysis we use the following values: 

TOfc, ^ 1.235GcV , fb, (2GcV) O.lSGeV^ . (38) 

where is obtained by using Eq. (A. 20) in Ref. 57] when assuming rrib^ = 1235 MeV. 

In order to solve all the five coupling constants gf_^^ and gf (1=1,2,3,4), we need to find five equations. Comparing 
Eqs. p4p and p7p . we find there are just five Lorentz structures which can be used to obtain five equations and solve 
these coupling constants. They are g^^, p'^Pu, Qf^Pv, p'pQv, and q^q,y Among them, p'^qi, and g^gj/ lead to the OPEs 
which are much larger than those obtained using others. Consequently, g^ and g^ are also much larger than other 
coupling constants. Using this fact we can simply our calculations a lot. 



1. Lorentz Structure 

First we choose the Lorentz Structure g^i,. After performing the Borcl transformation once, we obtain 

gS^ ^f^MjU^ f^-^ljT^_^^^^^^^^ (39) 

jg'JGGG) ^ 
384\/27r2 

We find that the OPE part only contains the condensate {glfGGG), and so the numerical result turns out to be quite 
small. Therefore, we will only keep term and omit other terms coming from gp'u'igpp' — ){guv' — ^"^i"' ), etc. 
to simply our calculation. The three-point correlation function p4p becomes 



r,^r^""Hp,p','z) = — '^t\f 2 2, X ( (40) 
g^..{gl. + g?{p ■ p') + .9?(p • q) + gfip' ■ q) + gfq^) 

D / / D / ^ D I D\ / / D I D\ I \ 

-5i PpPt, - .92 qi^^Pu - (.9i + .93 JPm'?" ^ (-92 + .94 )9m9!^ H ) ■ 

Using this three-point correlation function, we will show in the next subsection that g^ is zero because the OPE with 
the Lorentz structure p'^p',^ vanishes. Then after some calculations, we find that the coupling constant gf_^^ is around 
0.02 GeV, which is much smaller compared with the coupling constants g^ and g^ . 
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2. Lorentz Structures p^pt o-'^d, p'^j^qv 
For the Lorentz structures p'lJ)'^, and p'^q^, the OPE side vanishes, and so we simply obtain 

g? = 0,and53°=0. (41) 

3. Lorentz Structure q^iPi, 

In this subsection we choose the Lorentz structure Qf^p'^. After performing the Borel transformation twice, we obtain 



g,^Z}^,m^J,J,e . e . - ^^^^^ ^ [^^^ dT,\r,+ 



(42) 



glfGGG) /6riT2 + 4r| ^ d d ^ (4ri - 2T2)r^ 

"ri=l/T2,r2 = l/T| 



192V^7r2 V(ri+T2)3 ^ dr^ dr^' (n + ra)* 
In Fig.[3]we show the couphng constant obtained using this equation, which is around —1.2 GeV^^ 

3 I 5 1 3 

O 



I 



^0 



0.5 1 1.5 2 2.5 3 
Borel Mass^ [GeV^] 

FIG. 3: The coupling constant g2 as a function of M^, obtained using Eq. (|42p . The Lorentz structure is g/jpt- 



4- Lorentz Structure q^q^ 

In this subsection we choose the Lorentz structure g^Qiy. After performing the Borel transformation once, we obtain 

vy2 .^4; r7i2^-TO2^ V y 96^2^2 I28y27r2 ^ ' 

We can also perform the Borel transformation twice, and then obtain 

{g^ + g^)V2f.,mlj,J'^e-<J^h'<J^'^ (44) 

^ {glGG) / n 9 T1T2 \ (gf/GGG) / 6rf + 6t2t2 + 6tit2 

48V27r2V (Ti+ra)^ ^9ri 9r2 (n + T2)3 / 192V27r2 V {n + T2Y 



._d d , (-2t2 + 2rir2 - 2r2)TiT2 

'^dri dT2> {T1+T2Y 



ri = l/Ti^,r2 = l/T| 



The quantity g^ +54* is a function of the Borel mass, as shown in Fig. |4l The left and right figures are obtained 
using Eqs. (|43p and (H^ . respectively. We obtain that g^ + g^ is around 0.6 GeV~^, and so the coupling constant 
g2 itself is around 1.8 GeV~^. 
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BorelMass^ [GeV^] 



BorelMass^ [GeV"] 



FIG. 4: The quantity + as a function of M% . The left and right figures are obtained using Eqs. (|43p and (1 44 II , respectively. 
The Lorentz structure is g;jg,y. 



C. Decay Width 



The decay width of tti — )■ &i7r reads: 



2 X 



+ 2 X 



+ 2 X 



+ 2 X 



+ 2 X 



+ 2 X 



247r ml^ 
247r ml^ 
247r TO^rn^^ 
127r m^^ml^ 
127r mj^m^ 
127r mTTiW^^ 



(3 + ^) 



(2|g; 



^ + TO^^ + m 



\J\qA'^ +ml^^y\q^\'^ + ml + 



;^ |2 



'6i ' 



where gf'^^ is defined as 



= glrr + Vr<+W^92 + "4.9? = -0.58GeV . 



The decay width of tti — > biTT is around 3 MeV. 



V. THE DECAY MODE TTi ^ /i(1285)7r 
A. Three-Point Correlation Function 



(45) 
(46) 
(47) 
(48) 
(49) 
(50) 

(51) 



In this section we study the decay mode tti — > /i(1285)7r. The isoscalar axial- vector current for /i(1285) is 



and it couples to /i(1285) through [57 



f 1 - 



{0\jf-\hip,\))=mfjf,e'}^ 



(52) 



(53) 
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Then we can write down the three-point correlation function for the decay mode tti — >■ /i(1285)7r: 



X [g^^gp.'i'' + g?{{p ■ pOs^v - -P^'-P^O + 92 Hp ■ Q)9t^'^' - qt^'Pi^') 
+9f{{p' ■ q)9t^'u' - p[,'qu') + 9iiq'^9^c'ly' ~ qp.'qv')) ■ 

We note here that the current jjj^ can also couple to the pseudoscalar meson rj and 77' through 

{0\jl'Hp),v'ip))^iU,'P^- (55) 
Now the three-point correlation function is: 

T/,y Hp,P,q) - g^,p,(g^,ga+g^^P^)(gpa- ^) ^^^2 , p2)(^2 _ p,2)(^2 _ ^2) ■ (56) 

Only Eq. (154^ contains the Lorentz structures g^^, q^qi, and p'^^q^, which can be used to differentiate them. 
At the quark and gluon level, we obtain the following OPE which is divergent at — 0: 

i{gsqcrGq) ,3p ■ q p'-q. , i{qq){glGG) , p-q p-q p' ■ q 



ni^^---\p,p\q)=9,.^{P^^ 



6V2 P^q^ P''^q^ 18\/2 ^ P^q^ p'^q^ p"^q^ 



p'-q, ^ i{gsqcjGq){glGG) , 11,1 1 



'p'2g4^' 48^/2 ^ p'4g2 ' + (p/6^4 ) 3 (^^ ' 

( i{gsq<yGq) 1 , ^{qq){glGG) 1 ^ 

^'^^'^'^ V V2 pV 18V2 pV^ 

We note that the correction vanishes in this case. 

B. Numerical Analysis 

We use the following values to perform the numerical analysis (57l |: 

m/, = 1285MeV , //^ = 170MeV . (58) 

where //^ is obtained using the sum rules Eq. (4.52) in Ref. |53l when assuming m/^ = 1285 MeV. Comparing Eqs. (p4)) 
and p7p . we find several Lorentz structures, and we can obtain several different QCD sum rules. They are g^^^ p'^Pu, 
qfiP^, P'^qu, and qfj,q^. 

1. Lorentz Structure g^^i, 

First we choose the Lorentz structure g^jy. After performing the Borel transformation once, we obtain 



5*(5,gaGg)(.g,2GG) 1 



144V2 ^2 ■ 

Here the dominant condensate {gsqcrGq) vanishes after the Borel transformation, and so the coupling constant g^_^^ 
is calculated to be around 0.05 GeV, which is much smaller than g^ and g^ which we will study in the following 
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subsections. So we are in the same situation as the previous section. Similarly we omit some small terms, and the 
three-point correlation function (IMl) is simplified to be 



(60) 



^^'^ ' ^' (m2^ - p^){m\ ~ p'2)(m2 - q^) V 

5^496% + gfip ■ p') + 92 (p ■ q) + 9?ip' ■ q) + gfi^) 

D / / D / / D , D\ / / D , D\ \ I \ 

-ffi Pt,P^ - 92 q^iPu - (3i + 53 )P^Lq'y - [92 + 9a ht^qy) H j • 

2. Lorentz Structures p^pt and q^p'^, 
For the Lorentz structures p'^p'^, and qp,p'^, the OPE side vanishes, and so we simply obtain 

ffP = 0,and5°=0. (61) 

3. Lorentz Structures p'^^qv 

In this subsection we choose the Lorentz structure p'^^qv After performing the Borel transformation once, we obtain 

D ^Ui'^Uh'T^fif^ fe"™- 1 - e"'"'i ) = '^^^9sqoGq) 
^ - to2^ ' e 1 J 3\/2 ■ ^ ^ 

The coupling constant is a function of Mb, as shown in Fig. [5l The result is around 5 GeV"'^. 

20 rs ^20 

15 
10 




1 1.5 
Borel Mass- [GeV^] 

FIG. 5: The coupling constant g§ as a function of M%, obtained using Eq. (I62p . The Lorentz structure is p'^q^- 



4- Lorentz Structures q^qi, 

In this subsection we choose tlie Lorentz structure q^Qi^- After performing the Borel transformation once, we obtain 

jjV2f^,mlJf,mfJ' _2 /j^.j2 /j^,_f2 i{gsq<yGq) i{qq){glGG) 1 
ff4 2 5 ^ -e fi' ) = j= + —j= — . (63) 



V2 18\/2 T2 



The coupling constant is a function of Mb, as shown in Fig. [51 The result is around 7 GeV ^. 

C. Decay Width 

To calculate the decay width of tti — )■ /itt, we similarly introduce 



4'- = + (v/"^^, + l<?;iVm2 +|g;|2 + |^-;|2)^D ^ ^2^D ^ 2.41GeV. (64) 
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FIG. 6: The coupling constant as a function of M%, obtained using Eq. (|63p . The Lorentz structure is qy.qv 



The decay width reads 



247r toJ^ V m/j^ J 

^ (^''^^'^'' + "'^ + "^'/^ + 2yWT^vWT^) (66) 

fi 

S/ D 1^ |3 ^ 

y/i7ry4 I^ttI / /i-. 12 , _2 /Ir- 12 i ™2 i I;? |2 , _2 



127r mj^m^^ 



■(yWT^vW+^+i'r^i' + "^/,) (69) 



127r m'im'i 
^ II 



(2\q^\' + m^^ + m\ + 2y^i^pT^v^15rF+^) • (70) 



Here (|65|) mainly comes from the S-wave decay, (|66l) and (|67|) come from the D-wave decays, and the last three are 
the S wave and D wave interference terms. Numerically they are: 

dnH) = 23.2MeV , §§1 = 0.4MeV , §71) = O.SMeV , §^ = -2.6MeV , ^ = 3.6MeV , dZQl) = -1.2MeV , 

The decay width of vri — > /itt is around 24 MeV. 

Finally we want to note two important facts here. First, all the three Lorentz structures g^^, qfj,qi, and p'^j^q^ which 
lead to non-zero coupling constants are not contaminated by the other decay mode rin, as shown in Eq. (|56p . Secondly, 
from Eq. (|57p we find that all these three different Lorentz structures contain the condensate (ggqcrGq), the dominant 
one. However, for the Lorentz structure Q^^u: 

it vanishes after the Borel transformation. So the obtained gj^^ is much 
smaller than and obtained using the other two Lorentz structures, where this condensate does not vanish. 



VI. ISOSCALAR HYBRID 



In this section, we turn to study the decay modes of the isoscalar hybrid state ai with the quantum numbers 
jGjPC _ Now the 5-wave decay modes are ai(1260)7r and /i(1285)?7, etc., while the P-wave decay modes 

are 7777', etc. The interpolating current for ai is 

1 - A" 

= ^(""7.^' + d'^l.d'')^gsG;, , (71) 

and it couples to the isoscalar hybrid state CTi through 



(72) 
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The QCD sum rule leads to the following relations 

T^hT'^P,p\q) = r5^^^^-(p,p',q), (73) 

T^l.E^''^ip,p',q) = t{^cos9p-^smep)xT^l,-^f^^p,p',q), 
T^'pf^'(p,p',q) = *(^sin0p + y|cos0p) X T^^'^'^^{p,p',q) . 

From these relations, the decay width of ai — ^ ai(1260)7r, /i(1285)77, and -qr]' can be easily obtained. For ex- 
ample, with = 2.0 GeV and /^^ 0.013 GeV, lUa^ = 1.26 GeV, fa^ = 0.17 GeV, the decay width of 
CTi ai(1260)7r, /i(1285)77, and rji are around 770 MeV, 74 MeV and 0.3 MeV, respectively. 

VII. SUMMARY 

We have studied the decay properties of the hybrid states with J^'-^ ~ 1 ^ using the method of QCD sum rule. 
We have used the three-point correlation functions to extract the coupling constants, and then calculated the relevant 
decay widths. We work at the pion pole. First we take the tt mass in the denominator to be zero. Then we work 
at the limit — !■ and pick out the divergent terms only. Such a procedure simplifies the calculation greatly at the 
cost of throwing away all the finite pieces. 

TABLE I: The decay widths of the isovector hybrid state tti of l'^ J^^ — 1~1~^ and isoscalar one ai of l'^ J^^ = O'^'l"^ when 
the hybrid mass is taken to be 1.6 GeV, 1.8 GeV and 2.0 GeV respectively. 



Decay Modes 


Decay Angular Momentum 


Widths (MeV) 


M = 1.6 GeV M = 1.8 GeV M = 2.0 GeV 


TTl - 




P-wave 


180 MeV 


410 MeV 


800 MeV 


TTl - 


777r 


P-wave 


0.9 MeV 


2 MeV 


4 MeV 


TTl - 


T- rj'-K 


P-wave 


0.1 MeV 


0.4 MeV 


0.9 MeV 


TTl - 


^ blTV 


5+ -D- waves 


2.9 MeV 


14 MeV 


53 MeV 


TTl - 




5+ -D- waves 


24 MeV 


140 MeV 


410 MeV 


0"! - 


^ vv' 


P-wave 


< 0.1 MeV 


0.1 MeV 


0.3 McV 


CTl - 


T> aiTT 


S-t-D-waves 


60 MeV 


310 MeV 


770 MeV 


(Tl - 




S-f-D-waves 






74 MeV 



We have calculated the decay widths of both the isovector hybrid state tti with l'^ J^'-^ = 1^1 ^ and isoscalar one 
(Tl with I^J^'-' = 0+1 ^ . The present three-point correlation function formalism works well for the decay processes 
TTl — pn, TTl — rjir and tti — > tj'tt where the Lorentz structures are relatively simple. We have also discussed the 
modes biir, /itt. These two decay modes are complicated by the many Lorentz structures in the three-point correlation 
function, possible mixing between the S-wave and D-wave decay patterns, and possible contamination from other decay 
modes such as pir. We illustrate the variation of the decay width with the hybrid meson mass in Table HI The S-wave 
decay width of the bin, fin modes increases very quickly as the hybrid meson mass and decay momentum increase. 
But for the low mass hybrid meson around 1.6 GeV, the pn mode is one of the dominant decay modes. 

For the l""*" state 7ri(1600), the decay widths of ni — >■ pn, ni rjn and tti — >■ rj'n are around 180 MeV, 0.9 MeV 
and 0.1 MeV, respectively. The modes r]n,T]'n are strongly suppressed compared with the pn mode. Moreover, the 
pn mode is one of the dominant decay modes of the 7ri(1600), which is in strong contrast with predictions from some 
phenomenological models. 

In this paper we use the hybrid currents to study these hybrid states. We can also use the tetraquark currents 
which may also couple to these 1 ^ states. We have used such currents to study their masses, and found their possible 
decay modes through Fierz transformation [H, . We plan to study their decay properties like their decay widths 
etc in the future. This might be useful to know the internal structure of the 1 states. 

We suggest the experimental search of 7ri(1600) through the decay chains at BESIII: e^e^ -> J/ip{ip') — > tti -I- 7 
or e'^e~ — > J/ip{ip') ni + p where the tti state can be reconstructed through the decay modes ni ^ pn n^n^n^ 
or TTl — )• /i(1285)7r°. More details of possible experimental search at BESIII can be found in Ref. [70]. It is also 
interesting to look for ni using the available BELLE/BABAR data through the process e+e" — >■ 7* — pni, hini,^ni 
etc. Hopefully our result will be helpful to the experimental identification of the 1 ^ state. 
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Appendix A: The Decay Mode tti 6i(1235)7r With the Tensor Current 
1. Three-Point Correlation Function 

In this appendix, we consider the tensor current for the 6i(1235) meson and study the decay mode tti — &i(1235)7r: 

= d<Tf,^u , jjf = -^{u(Tf,^u - da^^d) , j™ = -Ocr^^d , (Al) 
which couples to 6i(1235) through 

(0|j™|&i(p, A)) ^ */™6^.„^6^/ . (A2) 
The three-point correlation function at the phenomenological side is: 



_ 10 ( _ Pf^th£_\( _ PaPa' \ 

^pcrafSP \9fJ.fJ.' 2 )\9aa' 2 / 



X {9bt-n9t^'y' +g?{{p-p')9f^'i^' - p't,'Pi^')) + 92 0? ■ q)9f^'i^' - qf,'P,y') 

+9f{{p' ■ q)9t^'y' -p'^'Qw) + 94iQ'^9t^'iy' - V?--')))) ■ 
Unfortunately, the tensor current can also couple to the vector meson p(770) through [6^ 



When the tensor current couples to p, the three-point correlation function is 



{0\jI^\p{P, A)) = ^f^{e'p. - etp,) . (A4) 



r™;™^^,^',?) = d^xd^ye^P'-e^''y{0\Tjjr{x)r\y)vUm) (A5) 



V nip nip / 

We note that the following relation exists between these Lorentz structures 

f-iiual3q°'p'^Pp + e^papq^P^Pf^ + f^ppapq°'p'^Piy = - ( ^lu paq" P^ + f^pupaP°'P ' 9 ■ (A6) 

And we can write Eq. (jA5P as 



rTiTE(PH)/ I \ TE a V^j/7ri ^Ti^ /p / A '7^ 

^ [ni^^~ p^){nij- p'^){nij^ - q^) 

The contribution of pi: appears in our analysis at the leading Lorentz structure ep^paap'", and makes it difficult to 
single out the contribution of bin. 
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At the quark and gluon level, we obtain the following OPE which is divergent a,t 



i{glfGGG) r 



, /a/' 2 i2\ 2r| \ ?,i{glfGGG) e^paaP" 



2 



\2. 



2. Numerical Analysis 

In the numerical analysis, we use the following value [6^: 

/b™(2GeV) = 180(20)MeV. 



(A9) 



The leading Lorentz structure e^paaP'" appears in the S'-wave decay. Both pn and &i7r modes contribute to 
this structure. If we naively assume that only the bin mode contributes, we obtain after performing the Borel 
transformation once 



mf 



■■"Til 

We can also perform the Borel transformation twice, and then obtain 

,,3 



7i{glfGGG) 
48V^7r2 



0^ ^/2f fTE./.-</T^g-m?,/T,^ _ z(.gg/GGG) -20r2 



ri = l/Tf,r2 = l/T| 



(AlO) 



(All) 



Using Eqs. (jAlOp and (jAllI) we can perform the numerical analysis. The extracted g^^^ is around 0.5 GeV as shown 
in Fig. [71 which is significantly larger than our previous result. We do not use this value to calculate the decay width 
since we are unable to subtract the contamination from the decay mode p-n very cleanly. One may use the extracted 
pTT coupling constant in the previous sections as input and subtract such contribution from the sum rule derived using 
the tensor current. However, such an approach is very crude with large uncertainty. 




1 1.5 2 2.5 
Borel Mass^ [GeV-] 



> 1.5 
O 

= 1 



' 0.5 



0^ 
0.5 



1 1.5 2 2.5 
Borel Mass- [GeV^] 



1 ^ 
1.5 

1 

0.5 





FIG. 7: The coupling constant gf_^^ as a function of M%. The left and right figures are obtained using Eqs. (IATo|| and (|ATT|| . 
respectively. The Lorentz structure is tppaap'" ■ 
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